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We investigate the spin chain model composed of periodic 
array of two kinds of spins Si and S2, which allows us to 
study the spin chains with impurities as well as the alternating 
spin chains in a unified fashion. By using the Lieb-Shultz- 
Mattis theorem, we first study the model rigorously, and then 
by mapping it to the non-linear sigma model, we extensively 
investigate low-energy properties with particular emphasis on 
the competition between the massive and massless phases. 

PACS: 75.10.Jm, 05.30.-d, 03.65.Sq 



Quantum spin chains have been providing a number of 
hot topics not only in condensed matter physics but also 
in statistical physics, quantum field theory, etc. In par- 
ticular, the competition between the massive and mass- 
less states in the spin chains has attracted much current 
interest. For example, doping impurities into massive 
systems such as the Haldane-gap system Q , the two- leg 
ladder system [^|J|] and the spin-Peierls system Q may 
cause the massless states, and sometimes stabilize the 
long-range magnetic order. For these systems it is an 
interesting issue to clarify how the quantum coherence, 
which produces spin gaps, is suppressed by impurities 
to result in a gapless state [Q. Also, if the impurities 
are magnetic and their concentration becomes high with 
periodic arrangement, the system naturally leads to the 
alternating spin chains which have also been studied in- 
tensively j^] . A common feature in these problems is how 
the massive and massless states compete with each other, 
providing a variety of interesting phenomena. 

In this paper we investigate the quantum spin chain 
composed of the periodic array of two kinds of spin Si 
and S2, putting particular emphasis on the formation of 
the massive and massless states. The model proposed 
here is related to the interesting topics mentioned above, 
and naturally interpolates the impurity models and the 
alternating spin chain models. For example, if we con- 
sider the case of dilute Si spins in the background of the 
host spins with integer 5*2, the model is related to the 
Haldane gap systems with magnetic impurities. On the 
other hand, in the high density limit of Si spins, it de- 
scribes the alternating spin chains. We will investigate 
low-energy properties of the model by using the Lieb, 
Schultz and Mattis (LSM) theorem and non-linear sigma 
model techniques. 



The model we will investigate consists of two kind of 
spins with nearest neighbor interaction. The Hamilto- 
nian is given by 
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(1) 



Sj, where 



with the periodic boundary condition Sj+n 
N (= MN') is the number of sites, which is assumed 
to be even integral in what follows. The spin quantum 
number Sj on each site is Sj — Si for j = 1 (mod M) and 
Sj = 5*2 for others: The periodic array of Si impurities 
is embedded in the host spins ^2 with the period M, i.e. 



Si®S2® S2 



S2 C^Si ® S2 



'S2. (2) 
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Therefore by tuning the period M, we can naturally in- 
terpolate the dilute impurity models (large M) and the 
alternating spin chains (Af being order of unity). The 
antiferromagnetic coupling Jj (> 0) is assumed to have 
a bond-dependence, Jj ~ J(l + 7) for j = and 1 (mod 
M; between spins with S = Si and S2) and Jj — J for 
others (between the same spins with S = 52). 

Let us start by specifying the ground state properties. 
Since spins are on a bipartite lattice, we can apply the 
Marshall theorem Q] to (|l|) : The ground state is specified 
by the spin quantum number S* = (IS*! — S2\N/M) for 
M — odd (even), which is non-degenerate except for the 
trivial spin degeneracy. Therefore the ground state of our 
model is either a spin singlet or ferrimagnetic. We are 
interested in the quantum effects on the spin liquid phase, 
so that we will mainly concentrate on the singlet cases 
(odd M) which possess a variety of interesting properties. 
We shall show the simple but remarkable fact that the 
universality class of the model is solely determined by the 
impurity spin 5*1, as far as the concentration of Si spins 
is finite. 

To begin with, we apply the LSM theorem [|j in or- 
der to address a question whether the above model can 
have a gapless excitation. Let be the M-sites trans- 
lation operator, which commutes with the Hamiltonian 
(|l|) by definition. Therefore, together with the Marshall 
theorem, we find that it acts on the ground state |$o) 
as r*^|$o) = e*'^|$o)- Next define the twist operator 
U = exp(2^ jS"!), and the corresponding twisted 
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state 1$) = U\^o)- The energy increment due to twist 
is easily calculated as {^\H\^/) - {^o\H\'^o) < const. /-^V- 
In order to ensure that the twisted state |$) is actually an 
excited state, we need to show the orthogonality of these 
states. We immediately find T^^C/T"*^ = (_i)2S„ff a/jj 
with 



S-eff = [Si + (M - 1)82] /M. 



(3) 



Note that the orthogonality condition is satisfied by the 
minus sign of the factor (— l^^SoffM^ Therefore, in the 
case 2SeeM = odd-integer, we can prove that the twisted 
state is orthogonal to the ground state ($o|'i') = 0, and 
hence the system has a gapless excitation, or alterna- 
tively, degenerate ground states. For example, this is the 
case for the spin chain with Si ~ 1/2 and ^2 = 1. In 
general, the existence of an excited state with 0(1 /N) 
energy depends on how many half-integer spins are in- 
cluded in the unit composed of M sites. 

Here we deal with the most interesting case, i.e. the 
integer S2 spin model with 5*1 spins doped, which may be 
related to the Haldane gap system with magnetic impu- 
rities. According to the above theorem, it is seen in this 
case: it is the impurity spin 5*1 that controls whether the 
system is massless (or massive with degenerate ground 
states) in the integer S2 background. For dilute impu- 
rities with half- integer Si, an excited state with energy 
of 0{1/N) ensured by the LSM theorem may be essen- 
tially the same as an impurity state in the Haldane gap 
systems. According to extensive work on the Haldane- 
gap systems, doping impurities induces the free S* = 1/2 
degrees of freedom Q near the edge of the valence-bond- 
solid states [ p^ , and they couple with the doped impurity 
spin 5*1, making a local object of the effective impurity 
state. If we increase the concentration of spins, 
correlations among these local impurity states become 
strong, and impurity bands are naturally formed by co- 
herent motion of effective impurity spins. Owing to the 
quantum effects, they again become massless (massive, 
see below) for half- integer (integer) Si , which will be dis- 
cussed below in terms of field theoretic methods. 

Based on the above rigorous results, we now construct 
low-energy effective field theory, which allows us to study 
the model with integer-5'i as well as half-integer-5i cases 
in a unified way. Moreover, we can see which possibility 
for the ground state is realized for half-integer-5'i cases, 
i.e. a massless state or a massive state with degenerate 
ground state. For this purpose, we shall map the system 
(0) to the non-linear sigma model [T^Jl^ by the use of 
the SU(2) coherent state path-integrals |14| . In what fol- 
lows, we again concentrate ourselves on the systems with 
singlet ground state, i.e, for the case M ~ odd. Note that 
the following analysis can be applied to the system with 
finite concentration of impurities in the thermodynamic 
limit. 

The coherent state in the spin-S* representation is 
here defined by |C) = (1 + ICl^)""^ exp(CS'-)15', m = S). 



Parametrizing C = tan^e^"^, we have (C|S^|C) = Sn, 
where n = (sin 9 cos cj), sin 9 sin (/>, cos 6) . By using the 
over-completeness relation / d/i(C)|C)(CI — 1 with the 
invariant measure (i/i(C) = "^^^^ , we can de- 

rive the path-integral representation of the partition 
function. Namely, by staggering the spin configura- 
tion as Sj = Sj{-~y^^n{j), the partition function Z = 
trexp(— /3i7) of the system (^ can be represented by 
Z = / IljLi ^/^[''^(i)] 6xp(— 5), where the action is ex- 
plicitly given by 
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s = zj2i-ysMMj)] 

-J2JjS,S,+i / dTn{j)-n{j + l). (4) 

Here ct;[n(j)] is the Berry phase acquired by the jth spin, 
= [ du f dTn{j) ■ [drn{j) x dun[j)] . (5) 

JQ Jo 

In what follows, we assume that the behavior of the field 
n(j) is similar to that for the uniform spin chain in the 
continuum limit. Therefore, we can divide n{j) into 
slowly varying part and fluctuation around it as usual, 
n{j) = m{j) + a{~y~^^l{j). This assumption should be 
confirmed to be valid by comparison with the results of 
the LSM theorem. 

Let us first calculate the Berry phase term in the con- 
tinuum limit, 

M N/{2M) 

^S=^E E [^l4,l-52(-)'=(l-4,l)] 
fe=l j = l 

X {uj[n{2Mj ~ M + k)] - Ljj[n(2Mj - 2M + k)]) 

•-^1 f f J2„ 



(6) 



with ScS defined in eq.(|^), where xi — x and X2 — t. 
To derive this formula, we have taken the continuum 
limit with respect to Ma lattice spacing: n{2Mj — M + 
k) - n{2Mj - 2M + fc) - Madim{2Mj - M + k) + 
{~Y2al{2Mj - M + k), where a is the lattice constant. 
This procedure implies that we are now concerned with 
the most important low-energy mode, although there are 
other massive excitation modes because the period of the 
lattice is M in our model. The validity of this procedure 
will be discussed by comparing the results with those of 
the LSM theorem. It should also be noted that the con- 
tinuum limit here corresponds to the limit N/(2M) 00. 
Namely, we assume in this approach the finite concentra- 
tion 1/M of the Si spins. Next, the interaction term is 
rewritten in the continuum limit, 



i— J J d xm ■ {dim x d2m) 
+iSca / / d^xl ■ (m x d2m) 
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M N/M 

Si = -JJ2Y1 [(1 + 7)51^2(4,1 + 4,m) 

k=l 3 = 1 

+Sl{l-5k,i){l-5k,M)\ 
X j dTn{M j - M + k)- n{Mj - M + k + 1) 



1 1 (fx [{dimf + Af] + const. 



(7) 



with (3 = Sl + 2[{l+-f)Si- S2]S2/M |T|. In this way, we 
end up with the action composed of the fields m{x) and 
l{x). The term of l{x) appears in a quadratic form, and 
can be easily integrated out, which consequently results 
in the following Lagrangian density 



C 



1 

2^ 



v{dimf + -{d2mf 



-—Cf^^m ■ {df_,m X d^m), 



(8) 



with 



9 = 2mSi (= 27riS'cffM), 
2 

S2 + {Si- S2)/M 
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S'cff 



V = 2 Ja 



Sl + 2[{l + -i)Si-S2]S2lM 
S2 + [Si - S2)/M 



(9) 



where S^s is defined in eq.(^. 

We recall here that the topological term with 9 can 
completely specify whether the system is massive or 
massless. By observing that only Si appears in 0, we 
again arrive at the conclusion that the universality class 
of the system is solely determined by the impurity spin 
Si, which completely fits in with the results of the ex- 
tended LSM theorem obtained above: for half-integer 
5i, the system is massless even if ^2 = integer because 
6 — iri (mod 27ri). Furthermore, for integer Si we can 
say beyond the LSM theorem that the system should be 
massive because 9 = (mod 27ri). By using the above 
formulae, we can discuss the interesting results deduced 
by the LSM theorem in more detail. The following state- 
ments are valid for the system with finite concentration 
of impurities, (a) Haldane-gap systems (5*2 = integer) 
become massless (are still massive) for half-integer (in- 
teger) Si impurities, (b) Massless spin chains (5*2 = 
half-integer) become massive (are still massless) for in- 
teger (half- integer) Si impurities. Especially in the lat- 
ter case, we would like to recall the analysis by Eggert 
and Affleck from which we can naively expect the 
following scenario: Integer Si impurities are screened by 
the two neighboring half-integer spins 5*1, forming local 
integer 5*2 — 25*1 objects. These effective spins may cou- 
ple with each other in a coherent way and produce a gap 
a la Haldane. 



We should mention here that our mapping to the sigma 
model may be justified for a high or intermediate concen- 
tration of impurities, but not for a dilute limit. However, 
a qualitative feature whether the system is massive or 
not is determined solely by the topological term, and is 
expected to be correctly specified via the present anal- 
ysis so far as the impurity concentration is finite in the 
thermodynamic limit. 

So far we have been mainly concerned with the impu- 
rity effects on spin chains. Here we discuss the opposite 
limiting case (high density of 5*1 spins), i.e. the alter- 
nating spin chain system in more detail. For example, 
when M = 2 with a half-integer 5*1 and an integer 5*2, 
our model ([|) neatly describes the alternating spin chain 
system which has a ferrimagnetic ground state jl^ , being 
consistent with experiments found so far . In this con- 
nection, modified alternating spin chains with a singlet 
ground state (quantum liquid phase) have been actively 
investigated, for which quantum fluctuations should play 
a vital role |18[l^. This problem may provide a new in- 
teresting paradigm of spin chains bridging the massive 
and massless Heisenberg chains. As a typical example 
JlOf , we here consider a slight extension of the Hamil- 
tonian (P with Sj = Si for j = 1 and 2 (mod 4) and 
Sj = ^2 for and 4 (mod 4), namely. 



5*1 (g) S'l (g) 52 5'2 5*1 5*1 



I 5*2 



(10) 



Antiferromagnetic couplings are Jj = J(l — 71) for = 1 
(mod 4; between S'l), = J(l — 72) for j = 3 (mod 4; 
between S2) and Jj = J[l + (71 + 72)/2] for j = others 
(between 5*1 and 6*2). According to the Marshall theo- 
rem, it is easily found that the ground state of this model 
is singlet. Though the system ( [l0| ) is invariant under four 
sites translation T'', the LSM theorem cannot apply for 
any values of Si and ^2, i.e, the twisted state is not or- 
thogonal to the ground state. This suggests that the sys- 
tem may be massive in general. However, we should also 
note that since the system has spin-alternation as well 
as bond-alternation, there could exist non-trivial mass- 
less phases between the massive phases, as is the case for 
uniform spin systems We will clarify these points 

using non-linear sigma model techniques. 

Keeping the above observations in mind, we derive the 
effective theory. We have the same Lagrangian (||), but 
with 

9 = 27riS'cff(l +7cff), 
2 



S'effx/l -7eff' 



V = 2Ja{S,s - AS^)Jl - 7^3, 



(11) 



where 



S,s = [Si + S2)/2, 

AS^ = {Si - S2)illSi - 7252)/(45eff), 
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7off 



{Si + S2Y - (5*1 - S'2)(7iS'i - 725'2) 



(12) 



We wish to study non-trivial cases, taking Si — 1/2 
and 5*2 = 1 as an example, and then generalize the 
discussions to arbitrary spin cases. To begin with, let 
us set 71 = 72 = 0. We find from eqs.(|ll|) and (|l^ ) 
that the topological term 9 is difi'erent from m (mod 
27ri), and therefore the system is massive. Even in this 
case, 7off ^ since it includes the effects of the spin- 
alternation. Let us next introduce the bond-alternation 
term 7, and observe what happens if the parameter 7 
changes from —1 to 1 for the case 71 = 72 = 7. At 
7 = — 1, the model ( |lO| ) becomes a set of isolated dimers, 
as seen from (p^). This system has massive excitations, 
which is indeed consistent with 9 = Q and g — > 00 in 
eq.(|l^). If we increase 7 up to 1, the effective 7eff changes 
from 7off = —1 to 1. The corresponding value of 9 
changes from to 6* = Zni. So it is predicted that we 
encounter a massless fixed point once aX 9 = i:i during 
this process. For the case 7 ^ 1, we need more careful 
treatment. As naively expected, in this case. Si and ^2 
spins strongly couple with each other antiferromagneti- 
cally and the model (^0|) should behave like the uniform 
5' = 1/2 spin chain with massless excitations (indeed we 
have 9 Sttz). We should note, however, that exactly at 
7 = 1 the system is separated into isolated pairs of spins, 
which may cause a singular behavior in the coupling con- 
stant, g ^ 00. Therefore, for 7—^1 the model exhibits 
a behavior quite similar to that for a half-integral spin 
chain although it is still massive. 

In general, for the case with half-odd integer Si and 
integer 5*2, we can predict that there appear massless 
phases 2>S'eff — 1/2 times as the bond alternation param- 
eter changes from 7 = —1 to 1. It is also shown that 
for the cases where Si and 52 are the same type of spins, 
i.e, half-odd integers or integers, there appear 2S'cff times 
massless fixed points. We wish to note that in the uni- 
form case 5*1 = 5*2 = 5eff our formulae reduce to those 
for the ordinary spin chain with bond-alternation, which 
has 2S'eff critical points for — 1 < 7 < 1 [Q. In this way, 
our approach based on non-linear sigma model techniques 
are quite powerful to systematically study the low-energy 
properties of alternating spin chain systems. 
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